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Abstract. The problem of planetary system formation and its 
subsequent character can only be addressed by studying the 
global evolution of solid material entrained in gaseous proto- 
planctary disks. We start to investigate this problem by consid- 
ering the space-time development of aerodynamic forces that 
cause solid particles to decouple from the gas. The aim of this 
work is to demonstrate that only the smallest particles are at- 
tached to the gas, or that the radial distribution of the solid 
matter has no momentary relation to the radial distribution of 
the gas. We present the illustrative example wherein a gaseous 
disk of 0.245 M 0 and angular momentum of 5.6 x 10 52 g cm 2 
s" 1 is allowed to evolve due to turbulent viscosity characterized 
by either a = 10 -2 or a = 10~ 3 . The motion of solid particles 
suspended in a viscously evolving gaseous disk is calculated nu- 
merically for particles of different sizes. In addition vve calculate 
the global evolution of single-sized, noncoagulating particles. 
We find that particles smaller than 0.1 cm move with the gas; 
larger particles have significant radial velocities relative to the 
gas. Particles larger then 0.1 cm but smaller than 10 3 cm have 
inward radial velocities much larger than the gas, whereas par- 
ticles larger than 10 4 cm have inward velocities much smaller 
than the gas. A significant difference in the form of the radial 
distribution of solids and the gas develops with time. It is the 
radial distribution of solids, rather than the gas, that determines 
the character of an emerging planetary system. 
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1. Introduction 

The combination of observational evidence with theoretical con- 
cepts strongly suggests that planetary systems, other than our 
solar system, ought to exist. We believe that planets form by the 
accumulation of solid matter 1 in accretion disks around young 

1 More precisely, the nonvolatile component of planets. Giant plan- 

ets, like Jupiter, formed as solid bodies to which a large amount of gas 

was subsequently added. 


low-mass stars. Thus, the seeds of the future planetary syst< n 
architecture are to be found in the radial distribution of solids 
Starling with this paper we embark on an investigation aiirm 
at calculating the time and space evolution of solid material ii 
protoplanetary disks. 

It is very important to realize that the distribution of sc it 
particles undergoes a global time evolution, which accomj a 
nies, but is not identical to, the global time evolution of th< 
gaseous component of the disk. The gaseous disk spreads an< 
its mass diminishes due to viscous torques. Initially, the solit 
particles are very small, so they are coupled to the gas, am 
both components evolve more or less concurrently. With thru 
the solid particles gain mass due to the process of coagulatum 
and decouple from the gas. The subsequent global evolution o 
solids differs from the evolution of the gas, the divergence bein: 
a function of coagulation process time scales and the magniti di 
of solid-gas interactions. Eventually, as the density of the ga- 
decreases and the mass of solid particles increases, the radial chs 
tribution of solids converges to its final form. This final forn i 
a major factor in determining the character of the nascent pi, n 
etary system, inasmuch as any further redistribution of m ;s 
is only due to the planet formation process via the growth c 
planetesimals into protoplanets (Safronov 1969; Greenberg e 
al. 1978; Kaula 1979; Welherill 1980; Nakagawa et al. 198 3 
Although this last process radically changes solid mass dis ri 
bution on the small-scale, it does not much alter the averc y 
large-scale radial distribution of solid material. 

To begin our discussion we elucidate the two essential c i 
cepts underlying our approach. First, we assume that prcio 
planetary disks are indeed accretion disks and that the accret e 
process is due to an anomalous viscosity. Astronomical obs:i 
vations of young low-mass stars have led to wide acceptanct c 
the idea that these stars are surrounded by accretion disks (fc 
recent reviews see Strom & Edwards 1993; Beckwith 19S4 
On the other hand, it seems that much of our solar system 
gross architecture has arisen from physical processes and c >r 
ditions characteristic for accretion disks (Levy 1993). Beca 
we believe that the gaseous solar nebula from which the ;c 
iar system evolved exemplifies contemporaneously observ: ni 
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disks around low-mass young stars, we can make a persuasive 
case that a// protoplanetary disks are in fact accretion disks. The 
conjecture is that some sort of anomalous viscosity operates in 
those disks. That is not necessary so, inasmuch as other pro- 
cesses such as centrifugal ly driven winds (Pudritz & Norman 
1983) may cause an accretion of mass as well. Nevertheless, 
we adopt here the “classical” viewpoint that accretion disks are 
powered by viscosity created by the interaction of turbulent ed- 
dies. We do not specify the source of this turbulence, but instead 
use the usual a parameterization introduced by Shakura & Sun- 
yaev (1973). 

Second, \vc accept the current planet formation concept ac- 
cording to which planets form by the accumulation of solid 
matter, beginning with dust, into progressively more massive 
particles, and eventually culminating in the few large bodies that 
comprise the planetary system (for a review see Weidenschilling 
1988). This is because, in the presence of turbulence, the for- 
mation of planets could not have occurred as a result of gravi- 
tational instability as was envisioned (Goldreich & Ward 1973) 
before the accretional (and by conjecture turbulent) character 
of protoplanetary disks was fully appreciated. In a hypothetical 
quiescent disk the metamorphosis of fine dust into roughly 1 - 
10 km sized solid planctcsimals occurs almost instantaneously 
once the dust layer becomes thin enough for gravitational insta- 
bility to take place. Thus the radial mass distribution of plan- 
etesimals directly reflects the radial distribution of the gas at 
the moment when the instability occurred. In a turbulent disk 
no such instantaneous aggregation of solids is possible. Instead, 
planetesimals are the final product of a relatively long-lasting 
process of solid material accumulation. Because the solid par- 
ticles are suspended in a permeating body of moving gas, they 
undergo radial transport while coagulating into progressively 
larger particles. Thus the radial mass distribution of planetesi- 
mals emerging from a turbulent disk does not reflect the radial 
distribution of gas at any particular epoch. In order to calculate 
the radial mass distribution of planetesimals, we have to follow 
the concomitant evolution of gaseous and solid components of 
a disk. 

Note that once both concepts presented above are accepted, 
and the coagulation time scale is not very short compared to the 
viscous time scale 2 , the above conclusion regarding the manner 
of calculating the mass distribution of planetesimals is unequiv- 
ocal. However, although both of our principal concepts are fa- 
vored by theory and observation, the notion that the radial mass 
distribution of planetesimals can be obtained directly from the 
surface density of gas (or vice versa as in so-called minimum 
mass solar nebula models) lingers in the literature. 

It is quite clear that an understanding of the global evolution 
of solid particles requires daunting calculations. Before entering 
into such calculations we need to see a way through some of the 
potential pitfalls. The first problem is the choice of the starting 
time and the choice of initial conditions. The processes leading 

2 Both coagulation and viscous time scales are neither global nor 
constant. Whal we need to compare here are local and momentary 
values of these time scales. 


to the accumulation of solids may, in principle, occur already 
in a parent molecular cloud even before a cloud collapses to 
form a star and a disk (Weidenschilling & Ruzmaikina 1994). 
However, it seems most reasonable to choose the onset of the 
disk’s viscous stage as the starting time for our considerations. 
The initial conditions for the subsequent disk evolution should 
be chosen to reflect conditions at the onset of the viscous stage. 
Because such conditions are poorly constrained, we would start 
from some arbitrary initial distribution of the surface density of 
the gas, restrained only by the values of the total mass, A/jj sk , and 
the total angular momentum, J disk . Fortunately, the specific form 
of initial conditions docs not much influence the evolution of the 
gas, inasmuch as the process is diffusive in nature and the initial 
conditions are forgotten after a time short in comparison with 
the evolutionary time scale. As for solid particles, we assume a 
certain initial size distribution. The initial space distribution of 
a solid matter follows the initial surface density of the gas. The 
subsequent evolution of solid material is not a viscous process, 
and there is a danger that the radial distribution of solids at later 
times is susceptible to the arbitrariness of the initial conditions. 
This problem is avoided if the solids are coupled to the gas 
during the time the gaseous component needs to forget the initial 
conditions. 

The second problem is the desired accuracy in describing 
the gas-solid interactions. Obviously, we have to consider aero- 
dynamic forces exerted on the solids by the gas because this is 
the only process leading to the radial evolution of disk’s solid 
material. However, we neglect all potential effects the presence 
of solids may have on the evolution of the gas. The major effect 
of solid particles on the evolution of the gas is through the opac- 
ity. In the protoplanetary disk the opacity is dominated by solids 
(dust grains). As the solids evolve due to coagulation and radial 
transport, so should the opacity (Morfill 1988). This effect is 
disregarded in our approach and the constant size distribution 
of dust is used as a parameter needed to calculate the opacity. 
Such an approach may, nevertheless, provide a reasonable ap- 
proximation to a disk’s opacity as long as there is some amount 
of small grains well mixed with the gas. Our reasoning is that 
the coagulation process is not 100% efficient, leaving behind 
enough small grains to support the opacity. The second effect 
of solid particles on the evolution of the gas is the gas-solid 
friction, which is usually negligible compared to other forces 
acting on the gas, and its omission is justifiable. 

The third problem arises because solid particles may either 
appear or disappear due to thermal processes. The temperature 
in the disk changes rapidly with the distance from the star. In 
addition, the entire disk cools down with time. The particles 
traveling inward may evaporate and new particles condense in 
regions that just cooled below a certain temperature. Particles 
with high velocities with respect to the gas may be obliterated 
by the heat generated by the friction even if the temperature 
of the ambient gas is low. The difficulty of keeping track of 
these processes is magnified by the compositional variety of 
solid particles, which may consist of “rock,” “metal,” or “ice,” 
as well as the combination of these basic elements. Because our 
intention is to concentrate on radial motion and coagulation as 
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the two fundamental processes governing the global evolution ot 
solids, wc neglect all above mentioned thermal processes. Thus, 
at least at this stage of our investigation, the solid particles are 
thermally indestructible, and no new particles can he created. 

This concludes the overall description of the approach wc 
have taken to calculate the global evolution of solids in pro- 
toplanetary disks. Perhaps, from such a model, we can read 
something about how sensitive the architecture of the resul- 
tant planetary system is to the initial conditions, such as 
and Jdiski as well as on assumed microphysics, such as a 
and coagulation slicking coefficient. The model relics on two 
dominant and interdependent physical processes, aerodynamic 
forces, which make the radial transport of solid particles possi- 
ble and coagulation, which makes the emergence of large bodies 
like planetesiinals possible. The interplay of time scales charac- 
terizing these processes produces the final radial mass distribu- 
tion of solids. In this initial paper wc focus on the aerodynamic 
forces. We have calculated space-time distributions of radial and 
transverse velocities of the various-sized solid particles relative 
to the accreting gas. We have also calculated the global evolu- 
tion of a single-sized, noncoagulating particles lor a number of 
particle sizes. The case of unisized particles illustrates the dif- 
ferences in time evolution between the gas and the solids. Using 
our results we are able to put some limits on coagulation time 
scales, based on the premise that the evolution of solid material 
must engender planets. 

In Sect. 2 we describe the evolution of the gaseous com- 
ponent of the disk. The methods used to calculate the effects 
of aerodynamical forces on solid particles are given in Sect. 3. 
In addition. Sect. 3 contains the results of specific calculations 
aimed at obtaining space-time distributions of radial and trans- 
verse velocities of solid particles. The methods and results of 
solid material evolution calculations are presented in Sect. 4. 
Finally, in Sect. 5 we summarize our findings and discuss some 
of the implications of our results. 

2. The evolution of the gas 

The evolution of the gas in a viscous protoplanetary disk is 

determined from the axisymmetric equations of continuity and 

motion. Cylindrical polar coordinates (r, </>, z) are used, with 
z — 0 corresponding to a disk midplane (equator). Because the 
disk is explicitely assumed to be turbulent, quantities such as 
gas density and velocity are broken up into mean and fluctuating 
parts: 

p = ~p + p and V = V + V ; 0) 

The equations are expanded out and the Reynolds averaging 
technique is applied to obtain equations describing the large- 
scale, mean behavior of the gas. The general forms of the mean 
equations for an axisymmetric disk with negligeable molecu- 
lar viscosity are given by Cuzzi et al. (1993). The effects of 
turbulence manifest themsel ves b y the presence of correlation 
terms, V/V', p'U/, and p'V-Vj in the large-scale equations. 
The presumpting that p* « p makes it possible to neglect 
the triple correlations p'V'V'. The gas velocity correlations are 


expressed in terms of turbulent viscosity v - (1 /3)aC s H , tur- 
bulent velocity - \fcxC Si and gradients ol thejnean velocity: 
vjv~l = V x 2 /?> and VJvJ = -y{d\ 7 l /dx j + dVj/dxi). In ex- 
pressing v and V x in terms ot the gas speed ot sound G s and disk 
half-thickness H we use the standard o-disk model and assume 
that the Rossby number is of the order ol unity. 1 he correlations 
of the form p’V' are modeled by the gradient diffusion hypoth- 
esis: /TV 7 = -DVp, with D being the coefficient of turbulent 
diffusivity. Canuto & Battaglia (1988) argued that D = i/Ko 3 * * * . 
where Ko = 1.2 - 1.7 is the Kolmogorov constant. 

As wc arc concerned here with the evolution of the radiai 
distribution of the gas, wc assume V z = 0 and ignore the ver 
tical dependence of physical quantities describing the state oi 
the disk. Under such assumptions the only nonvanishing corre 
lations are: 



The r component of the equation of motion determines 
large-scale tangential velocity. We anticipate that a protoplan 
etary disk is geometrically thin. In_a thin^accretion disk th< 
large-scale velocities are such that V# » V r ~ vjr and the 
component of the equation of motion reduces to: 

V = 

r r p dr p dr \ 3 J 

_2Dd^d^_K±( rD dp\ (5 . 

~p dr dr ~pr dr \ dr ) 

Here V v = ( GM/r ) l/2 is the Keplerian velocity and M is th 
mass of the central star. Of the five terms on the right-hand sid 
of Eq. (5) the first, due to centrifugal acceleration, is the larges 
The second, due to the gas pressure gradient, is smaller hy 
factor of the order of ( H/r ) 2 . The third, “turbulent pressure ’ 
term, is smaller by a factor of the order of (V K /C S ) (H/r) 
The two last terms, due to turbulent_ diffusion, are smaller b 
a factor of the order of a(H/r) 2 (\^/V k ). Thus the existence 
of sub- sonic turbulence does not influence the mean tangenti I 
velocity of the gas, which is slightly smaller than the Keplerian 
velocity because dP/dr < 0. This small difference is crucial 
for the gas-solid interaction, but unimportant for the evolution 
of the gas, where V# = V k approximation is sufficient. 

The mean radial velocity of the gas is computed from the n 
component of the equation of motion. At the lowest order, th s 
component can be written as 

The convective mass flux based on an average radial veloci y 
(the first term on the left-hand side of Eq. (6)) is, in general, of t! e 
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same order of magnitude as the diffusive mass 11 ux (the second 
term on the left-hand side of Eq. (6).) Together, they transport 
as much mass as the conservation of angular momentum allows 
(the right-hand side of Eq. (6)). 

The Reynolds averaged continuity equation in the absence 
of sources is 


dp 

0t 


( rpV r + rp'V;) 


r dr \ 


pV r - rD- 


dp 


dr 


(7) 


The evolution of a gaseous disk proceeds according to well 
established principles: the surface density decreases as the mass 
is lost to the central star (sec Figs. 4 and 5), the temperature 
decreases as well, and the disk spreads as a certain portion of 
the disk’s mass moves outward, carrying angular momentum. 
The physical quantities describing the slate of the gas change 
dramatically during the evolution, thus constantly changing the 
aerodynamic regime encountered by solid particles suspended 
in the gas. 

3. Aerodynamic forces on solid particles 


Substituting V r from Eq. (6) into Eq. (7) and integrating over the 
2 coordinate we obtain the familiar equation for time evolution 
of the surface density (E = f pdz) of the gas 


3E _ 3 d ( 
dt r dr \ 


1/2 if 

dr 




M 1 - 


( 8 ) 


Note that Eq. (8) does not depend on turbulent diffusivity, as it 
only encapsulates conservation of mass and angular momentum. 
Because turbulent viscosity, u , is not an explicit function of time, 
but instead depends only on the local conditions within a disk, v 
can be expressed as v - jv( E } r) and Eq. (8) can be solved subject 
to boundary conditions on the inner and outer edges of a disk 
and the opacity law. Given E(r, t) } we can algebraically find all 
other disk variables. In particular, V r can be calculated from Eq. 
(6). Beware that the average radial velocity based on Reynolds 
averaging does not give an accretion rate, or — 2?rrE V r / M. In- 
stead, the sum of convective and diffusive fluxes give the proper 
accretion rate, M = -27rrEV^ + 2 tt rD(9E/3r). 

The methods of how to solve Eq. (8) and how to calculate 
other physical quantities describing the state of the protoplan- 
ctary disk are given in papers by Ruden & Pollack (1991) and 
Reyes-Ruiz & Stepinski ( 1 995). Here we follow the methodol- 
ogy of Reyes-Ruiz & Stepinski 3 . In our numerical calculations, 
we start from the initial surface density profile of the form 


E(£ = 0,r) = E 0 [l +(r/r 0 ) 2 ] 5 (9) 

This distribution corresponds to E being practically uniform 
between the inner radius, which we have chosen to be r m = 
0.037 AU, and ro- The particular choice of r tn is not important 
as long as it is much smaller than the size of the disk. Values of 
ro, Eo, and s arc calculated to correspond to the desired initial 
disk mass and its angular momentum. For calculations presented 
in this paper we started with a disk of 0.245 Af© and an angular 
momentum of 5.6 x 10 52 g cm 2 s -1 . These correspond to an 
initial surface density distribution characterized by Eq = 10 4 
g cm' 2 and ro = 15 AU, initial conditions that correspond 
closely to those considered “standard” in protoplanetary disk 
evolution calculations (Ruden & Pollack 1991). We examine 
disk evolution for two values of a: a = 0.01 and a = 0.001. 


The aerodynamics of solid bodies in the context of the solar 
nebula was first studied by Weidenschilling (1977). However, 
whereas his work pertained to a quiescent gas environment, we 
consider the aerodynamics of solids particles in a turbulent disk. 
The presence of turbulence complicates the problem because of: 
a) the existence of the mean radial velocity of the gas, b) the 
diffusion of solids due to collective action of turbulent eddies 
advecting particles in all directions. In accordance with the prin- 
ciples outlined in Sect. 1, we not only consider a turbulent disk, 
but we calculate the space-time variations of gas conditions in 
a self-consistent manner. We want to investigate the global evo- 
lution of solid bodies with sizes from ~ pm to few hundreds of 
meters embedded in the accreting turbulent gas. The basic time 
evolution equation for a “fluid” of solid particles of mass in is 
given by 


d\ d 

Pd-fa = -Pd^y ~ Pd 


= -PdV'I'-pd 


Vrf- V 

t s (m) 


t s (m) 


(10) 


where $ is the gravitational potential, v is the relative velocity 
between a particle and the gas, and t s is the so-called stopping 

time defined as t s (m) = mv/F D . Here v = |v| = yjv* + u 2 , and 

F d is the drag force, which depends on the size of the particle, 
the density of the gas, and the speed of the particle relative to the 
gas. Subscript d denotes quantities describing the state of the 
solids 4 , whereas quantities describing the gas have no subscript. 
We now proceed in a fashion similar to that described in Sect. 2 
for the gas. Particles density and velocity are divided into mean 
and fluctuating parts: 


Pd = Pd + Pd and V d = \ d + \' d 


( 11 ) 


Radial and tangential components of Eq. (10) are expanded 
out as we seek to calculate V r d and V^. We follow Cuzzi et 
al. (1993) in expressing correlation terms as V^V- /Sc, 

where Sc is called the Schmidt number and is given by the ex- 
pression 


Sc = (1 +O k i s ) 



( 12 ) 


3 Reyes-Ruiz & Stepinski include the magnetic field in their con- 
siderations. Here we are using their numerical code but we explicitly 
exclude magnetic fields. 


4 In this paper solid particles are called particles, solids, dust, or 
grains interchangeably. We have chosen subscript d for “dust” to denote 
all quantities describing solid particles. 
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Table 1. Summary of drag laws used in our calculations. 


Reynolds number 

Stokes law 

Epslcin law 

Re = Ispvfi) 

A < 4s/9 

A > 4.s/9 

Re < 1 

/•'n = 24/ic - 1 /•' 

Fd = 4 -fp* 2 iC\ 

1 < lie < 800 

F n = 24 Re-^F 


Re > 800 

F 0 = 0.44 F 



f=™-,4 



Here Q k denotes the Keplerian angular velocity, and the turnover 
lime of turbulence is approximately given by 1 /Qy as we have 
assumed that the Rossby number is of the order of unity. Parti- 
cles are imperfectly coupled to turbulent eddies because: a) the 
lifetime of an eddy can be shorter than the slopping time of a 
particle (this is described by the first factor on the right-hand 
side of Eq. (12)), b) a particle has a systematic relative velocity 
with respect to the gas and can cross an eddy in time shorter 
than the slopping time (this is described by the seco nd fac tor 
on the right-hand side of Eq. (12)). The correlations are 
expressed as -D j and the diffusion coefficient is given by 
Dj = D/Sc (see Cuzzi et al. 1993, as well as Dubrulle et al. 

1995). 

Taking advantage of the fact that V rd « C s we can write 
the r component of Reynolds averaged Eq. ( 10) as 

%d = Yl + J_ JL ( E¥l. \ + TL 

r r ~pd dr \ 3Sc / t s + 

2D (13) 
pdSc dr dr ’p d r dr \ Sc dr ) 

When averaging the gas-solids coupling term we assume that 
p d y/t s ~ ~pd v/£**, where t sm = t s (v) if v > V t and t s * - t s (V t ) 
if u < K. In addition, we have assumed that p^V/^-p^V/ ^ 0. 
Of the five terms on the right-hand side of Eq. (13), the first 
is the largest. The second, fourth, and fifth terms, all due to 
coupling of particles to turbulent eddies, are negligeable for 
subsonic turbulence. Thus, V^ d differs slightly from Vy because 
of Tv/£ s * term. 

Averaging the ^-component of Eq. ( 10) we can ca lcula te the 
mean radial velocity of particles. Using the fact that V$ d ~ Vy 
wc obtain at the lowest order 


if 


2 p d 
£lyt s + 


(14) 


Note the formal similarity between Eqs. ( 1 3)-( 1 4), describing 
the motion of particles, and Eqs. (5)-(6), describing the motion 
of the gas. Due to the coupling of particles to turbulent eddies 
of the gas, they experience “pressure” (the second term on the 
right-hand side of Eq. (13)), as well as “viscosity” (the first 
term on the right-hand side of Eq. (14)). In addition there are 
differences between mean velocities of particles and the gas 
resulting in the drag force (terms proportional to v in Eqs. (13)- 

(14)). 


3. 1. Drag laws 

To close the system of Eqs. (13) and (14) we need to spcci > 
t s (m, v, gas) = mv/F 0 . We assume that a solid particle is a 
sphere of radius 5 and bulk density p s> thus s = (2>m/4np s ) ] ' 
The drag regime depends on die size of the mean free pa f 
of gas molecules, A, as compared with s. Table 1 shows t c 
forms of different drag laws used in our calculations. The g is 
molecular viscosity rj is given by the formula p - p\C s /2, c 
the Reynolds number can be expressed as Re = 4(u/C s )(.s/ , ) 
As wc do not expect that particles have supersonic velocit; :s 
with respect to the gas, the Epstein regime is tantamount u 
Re < 1. In the Re < 1 regime the transition from Stok ^ 
law to Epstein law occurs for X/s = 4/9. We calculated i.n 
mean free path of gas molecules from the usual formula A - 
I /ncro, where n is the number density of gas molecules and . 
is the total scattering cross section. The major constituent r 
die gas in protoplaneiary disks is molecular hydrogen. Taki =g 
the geometrical approximation to the cross section we obi; i 
A = ??ih 2 /7t/97’ } t. Because the state of the gas is known (s ■( 
Sect. 2) the mean free path can be calculated at any time and t o; 
any disk location. This information combined with the data u 
particle mass and relative speed (or V x if v < V t ) sets the vai it 
of the Reynolds number and the drag regime, thus allowing i lit 
closure of the system of Eqs. (13) and (14), albeit in the impli *i 
form, as knowledge of a particle’s relative speed is requii *,c: 
to establish a drag regime. This implicitness of our equatiem 
is not a serious problem inasmuch as the solution is obtain ;c 
numerically anyhow and an appropriate algorithm is applied a 
achieve the self-consistency of the solution. 

3.2. Space-time distribution of solid particles mean velocitit > 

Eqs. (13) and (14) constitute the closed system for two i n 
knowns, Tv and because V$ d can be eliminated using <h< 
relation V^~ d = 0$ + % with % - V k + (r/2pV k )(dp/dr) am 
V rd can be eliminated using the relation V rd = Tv + V r with K 
given by Eq. (6). With such an elimination, and keeping o: h 
terms linear in Tv and TJ^, we obtain 



v r 

1 

dP 1 


2vj, 


P 

dr n k 


and 





W + 

2v^ 

-V T 

— 3r~ i r 




dr V Sc / 


(D - 2u) 1 dpi 
Sc p2 dr 


The time evolution of Tv and T^ can be found for particles c 
any given size by solving the system of Eqs. (15-16) as al 
coefficients in these equations depend on properties of the ga 
and thus are known functions of space and time. The last tern 
in Eq. (16) is the exception as it depends on the density o 
the particles. The space-time evolution of particles density i 
discussed in Sect. 4, here we remark that calculating parti 1 
densities requires the knowledge of Tv and v$. Therefore i 
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general, Eqs. (15-16) have to be solved simultaneously with 
(he particles density evolution Eq. (18). Note, however, that the 
last term in Eq. (16) vanishes if D = 3m, which is consistent 
with estimation of D by Canuto & Battaglia (1988) (see the 
begining of Sect. 2). In the bulk of our calculations presented 
in this paper we assume D = 3m. This permits the discussion 
of particles mean velocities independently from the discussion 
of their densities. We have also performed calculations for D = 
is, which requires the solution of a coupled velocity-density 
problem, and have found only expected, quantitative differences 
between the results for D = v and the results for D = 3m. Thus, 
it seems that although the condition D = 3m is computationally 
advantageous, there is nothing physically unique about it. 




Fig. 1. Radial distributions of V Td , the mean radial velocity of a solid 
particles; V r% the mean radial velocity of the gas; and v^ t the mean 
relative transverse velocity between particles and the gas, respectively. 
The times shown are t = 10 4 yr, t = 10 5 yr,and* = 1 0 6 yr. Particle have 
all the same radius s = 0. 1 cm. Solid lines indicate a positive velocity 
(an outward movement for radial velocities and particles faster then 
the gas for uj), whereas dash-dotted lines indicate a negative velocity. 


J-ct’s consider the gas evolution scenario described in Sect. 
2, which starts at t — 0 with a disk ol 0.245 jV/^ and an angular 
momentum of 5.6 x 10 52 g cm 2 s“ ! . The evolution proceeds 
with dimensionless coefficient of viscosity a ~ 0.01. Our goal 
is to calculate the time evolution of V rd and V d)d for particles 
of different sizes embedded in an evolving gaseous disk. The 
knowledge of particle mean velocities is important inasmuch as 
these velocities capture the bulk of particle dynamics. Beware, 
however, that the contribution of turbulent diffusion to particles 
dynamics is crucial for explanation of certain features of their 
evolution. First, we have found out that very small particles, 
with 5 < 0.1 cm, can be considered perfectly coupled to the 
gas. Their mean velocities match the mean velocity field of the 
gas. On the other hand, very large particles, with s > 10 5 cm, 
can be considered completely decoupled from the gas. Their 
mean velocities remain practically unchanged with lime. Parti- 
cles with sizes between these two extremes are partially coupled 
to the gas. 

1 stl0ws time evolution of V r (mean radial velocity of 
the gas), V rd (mean radial velocity of a panicle), and (tan- 
gential component of mean relative velocity between a particle 
and the gas) for particles with s = 0. 1 cm. For the purpose of 
the present discussion we label such particles “small.” The tan- 
gential velocity of the gas is not shown in Fig. 1 as it always 
stays very close to theJCeplerian velocity. As we have already 
mentioned in Sect. 2, V r is not always the accurate indicator of 
the mass flux. The distinctive feature of the mass flux distribu- 
tion in the gaseous accretion disk is the existence of a stagnation 
radius that separates the region of gas accretion from the region 
of gas decretion. This stagnation radius moves outward during 
disk evolution. Examining Fig. 1 we notice that although V~ r has 
the same sense as the mass flux throughout the entire accretion 
zone and in the inner decretion zone, it is oriented inward in 
the outer decretion zone even so the mass flows outward there. 
This is because the transport of mass in the outermost disk, 
where the steep gradient of density is maintained, is dominated 
by turbulent diffusion. Inward convective mass flux is neces- 
sary to reduce the outward diffusive mass flux so that Eq. (6) is 
satisfied. 

Small particles remain rather strongly coupled to the mean 
flow of the gas; t5J is very small and V, ~d » K- Throughout the 
accreting part of the disk the magnitude of vj remains within 
the same order of magnitude, about 10" 3 cm/s at t = 10 4 yr, 
10~ 2 cm/s at t = 10 5 yr, and finally KT 1 cm/s at t = 10 6 yr. In 
comparison, the Keplerian velocity in this region of the disk is in 
the range of I0 6 — 10 7 cm/s. The particles have super-gas (faster 
than the gas) mean tangential velocities and thus experience, on 
average, a tangential head wind. The particles have small super- 
gas mean radial velocities {V r d < V r < 0) and thus experience, 
on average, a radial head wind as well. Throughout the outer, 
decreting part of the disk, the particles are less coupled to the 
mean flow of the gas. The magnitude of v$ changes by orders 
of magnitude and can be as large as 10 4 cm/s at the outer edge 
of the disk. For times up to few x 10 4 yr, the particles in the 
inner part of the gas decretion zone move on average outward 
with slightly sub-gas velocities, thus experiencing a radial tail 
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Fig. 2. Radial distributions of V r d, V r , and v#, respectively. The times 
shown are t = 2 x 10 4 yr, t = 6 x 10 5 yr, and t = 10 5 yr. Particles have 
a radius s = 10 cm. Solid lines indicate a positive velocity (an outward 
movement for radial velocities and particles faster then the gas for t?^), 
whereas dash-dotted lines indicate a negative velocity. 

wind. However, in the outer part of the gas decretion zone, the 
average radial velocity of small particles is directed inward, as is 
the average radial velocity of the gas. Nevertheless, mass fluxes 
of both particles and the gas are outward, as they are dominated 
by turbulent transport. 

Fig. 2 shows time evolution of V r > V r di and ^ f° r particles 
with s = 10 cm. We label such particles “medium.'’ Medium 
particles are much less coupled to the mean flow of the gas than 
small particles: v$ for medium particles is orders of magnitude 
larger than for small particles. There is not much change 
in the magnitude of vj throughout the disk evolution. Medium 
particles move on average with super-gas tangential velocities 
everywhere in a disk and any time during the disk evolution. 
Also, medium particles always move on average inward. Note 
that |Kd| >> |V^|, thus these particles are “fast” in their ra- 


Fig. 3. Radial distributions of V r d> Vr> and respectively. The time 
shown are t = 10 4 yr, t = 10 5 yr, and t = 10 6 yr. Particles have • 
radius s = 10 4 cm. Solid lines indicate a positive velocity (an outwar i 
movement for radial velocities and a particle faster than the gas for u £ 
whereas dash-dotted lines indicate a negative velocity. 


dial movement. With average radial velocities of the order c f 
10 3 -10 4 cm/s, medium particles have evolutionary time scale 
comparable to 10 4 yr. Therefore, handling of the medium part 
cles encounters the problem described in Sect. 1. As their eve 
lutionary time scale is comparable to the time the gaseous dis . 
needs to forget the arbitrary initial conditions, their propertit 
at later times may be sensitive to the initial conditions. To cii 
cumvent this problem we introduce medium particles into th ■ 
gas at t = 10 4 yr, instead of t = 0 as we have done for the sma 1 
particles. Also, we finish our calculations for medium partick s 
at t - 10 5 yr, as medium particles are severely depleted afti r 
this time. In fact, there would be no medium particles left in 
the disk after this time if not for the contribution of turbulent 
diffusion to the dynamics of particles. 
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Fig. 3 shows time evolution of V r , V riU and for particles 
with s = 1 0 4 cm. We label such particles “large ” Large particles 
arc only weakly coupled to the mean gas ilow: uj remains rela- 
tively large as particles move with almost a Kcplcrian velocity, 
faster than the gas that is pressure supported. Large particles, 
much like the medium particles, move with super-gas tangen- 
tial velocities everywhere in a disk and at any time during the 
disk evolution. They also always move inward. Note, however, 
that conversely to the case of the medium particles, for t > 10 5 
yr, \V r d\ « IK-1, so the radial movement of large particles is 
insignificant, and becomes more insignificant with time. Thus, 
large particles arc “slow’' in their radial movement. They remain 
on Keplcrian orbits drifting very slowly inward. 


4. Global evolution of single-sized, noncoagulating particles 


In order to visualize the global effects of solid particles decou- 
pling from the gas, it is convenient to calculate the evolution 
of single-sized, noncoagulating particles. This is a hypotheti- 
cal evolutionary scenario, by itself not likely to be pertinent to 
establishing the architecture of an emerging planetary system. 
Nevertheless, it illustrates the difference in time evolution be- 
tween the gas and the solids. It provides a valuable point of 
reference for more realistic calculations. In the case of D ^ 3u, 
computation of space-time distribution of ~pd is necessary to es- 
tablish a particle velocity field (see Sect. 3.2.) Here, however, 
we assume that D = 3v and velocity of particles are known 
beforehand. 

The ensemble of identical particles can be treated within the 
framework of hydrodynamics, or as a “fluid” with “pressure” 
and “viscosity” provided by coupling between particles and tur- 
bulent eddies of the gas. By analogy with the gas, the evolution 
equation can be obtained from mass conservation supplemented 
by the <p component of the equation of motion. The Reynolds 
averaged continuity equation for fluid of single-sized particles 
in the absence of sources is 


m 

at 


~l^{ rpdVrd+rp d V rd) 

1 9 ( —JJ- D dpj\ 


(17) 


Substituting Ka from Eq. (14) into Eq. (17) and integrating over 
the z coordinate we obtain the equation for time evolution of 
Ed, the surface density of solid particles 



One can also define Md as the accretion rate of a solid material 
through a disk. As in the case of the gas, the mass flux of solid 
particles is the sum of convective and diffusive fluxes. Unlike 
Eq. (8), which directs the evolution of the gas, Eq. (18) is not 
a diffusive-type equation. Instead, it is an advection-diffusion 
equation. Because Eqs. (8) and (18) are of different types, a 


single numerical method that would concurrently advance both 
of them in time is difficult to contrive. Taking advantage of 
our assumption that the evolution of the gas is not affected by 
the presence of solid particles, we first solve Eq. (8) using an 
implicit scheme, then obtain the space-time distribution of v# 
by methods described in Sect. 3, and finally solve Eq. (18) using 
the operator splitting method. In such a method the advcclive 
term in (18) is treated by numerical method of characteristics, 
whereas an implicit scheme is applied to the diffusion term. 

In general, particles can be divided into “long-lived” (LL), 
characterized by evolutionary time scales comparable to or 
longer than the evolutionary time scale of the gas, and “short 
lived” (SL), characterized by time scales shorter then the life- 
time of the gaseous disk. The small particles, which are strongly 
coupled to the gas and thus evolve on an approximately viscous 
time scale, and large particles, which are decoupled from the 
gas and move toward the central star very slowly, are LL parti- 
cles. The medium particles have radial velocities larger than the 
radial velocity of the gas, they accrete quickly onto the central 
star, and therefore are SL particles. The terms LL and SL refer 
to an average lifetime of a particle. Because of the existence 
of a random component in particle motion, few individual SL 
particles may survive for a long time, and some individual LL 
particles may have a short lifetime. 

Fig. 4 shows an evolution of the surface density for LL 
particles. Our computation starts at t = 0 with a disk of 
0.245 Mq and an angular momentum of 5.6 x 10 52 g cm 2 
s" 1 . Evolutionary scenarios with a = 0.01 and a = 0.001 
are calculated. The initial surface density of solid material is 
Ed(r, t = 0) = 10~ 2 E(r t t = 0) to reflect that about 1% of disk 
mass is in the form of solid particles. We follow the evolution 
of LL particles up to t = 10 6 yr. Solids with a radius of s = 0. 1 
cm represent small particles, and solids with a radius of 10 4 cm 
represent large particles. 

If the entire mass of the solid material is concentrated in 
particles smaller than 0.1 cm, then solids are perfectly coupled 
to the gas and the radial distribution of solids follows the radial 
distribution of the gas as indicated on Fig. 4 by dotted lines. For 
the purpose of the present discussion we refer to such a distri- 
bution of solids as the putative distribution (it differs from the 
distribution of gas by a factor of 100). Now assume, however, 
that the solid material is concentrated into 0.1 cm particles. It 
is clear from Fig. 4 that such solids slowly decouple from the 
gas. First, let’s discuss the a = 0.01 scenario. At t = 10 4 yr 
the decoupling of small particles from the gas is not yet pro- 
nounced. The decoupling becomes noticeable at £ = 10 5 yr, 
when solids located at the outer edge of the disk spread outward 
at the slower rate than the gas. This, combined with the fact that 
the inward radial velocity of solids is slightly larger than the 
radial velocity of the gas (see Fig. 1), causes the surface density 
of solids within the inner 10 AU to be enhanced with respect to 
the putative distribution. At t = 10 6 yr the difference between 
the surface density of small particles and putative surface den- 
sity increases. Within the inner ~ 10 AU the mass density of 
particles is enhanced with respect to the putative density, but 
ouward of ~ 10 AU the mass density of particles is actually 
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Fig. 4. Radial distributions of surface density of solids at selected times (£ = 10 4 yr, t = 10 5 yr, and t = 10 6 yr) for LL particles. The left-ha .c 
panels correspond to the evolutionary scenario with a = 0.01 , whereas the right-hand panels correspond to the scenario with a = 0.001 . Dott ;d 
tines show the surface density of the gas x 1(T 2 . Solid lines show the surface density of 0.1 cm particles and dashed lines show the surface 
density of 10 4 cm particles. 


depleted as compared to the putative density. The a - 0.001 
shows a similar behavior of small particles. 

Now assume that the solid material is concentrated into 10 4 
cm particles. Examination of Fig. 4 reveals that such large par- 
ticles are largerly decoupled from the gas. In the a = 0.01 sce- 
nario the radial distribution of surface density of such particles 
in the inner disk shows relatively little change during t = 10 6 
yr of evolution. Therefore, a large difference between the mass 
density of large particles and the putative mass density develops 
with time in the inner portion of the disk. In the outer portion 
of the disk large particles evolve faster than the putative mass 
density, again resulting in build up of a significant difference 
between both distributions with time. Overall the radial distri- 
bution of large particles steepens with time, whereas the putative 
mass distribution flattens. In the a = 0.001 scenario, after some 
initial time, large particles evolve faster than the putative mass 


distribution. This is because the a = 0.001 disk evolves me x 
slowly than the a = 0.01 disk and its gas density remains hi J; 
for much longer. That, in turn, makes the drag force on large p ar- 
ticles stronger and longer lasting, thus enabling them to evol ve 
faster. Again, a large difference between the mass density oi 
large particles and the putative mass density develops with tin e 

Fig. 5 shows an evolution of the surface density for SL p j 
tides. Our gas computation starts at t = 0 with a disk of 0.2 +‘ 
Mq and an angular momentum of 5.6 x 10 52 g cm 2 s“ 1 . Evo u 
tionary scenarios with a = 0.01 and a ~ 0.001 are calculate d 
To eliminate the influence of arbitrary initial conditions on tnt 
evolution of SL particles, for which evolutionary time sea e 
are comparable to 10 4 yr, we start our particle computation a 
t = 10 4 yr with E d (r,t = 10 4 ) = 10- 2 £(r,f = 10 4 ) and v, 
follow their evolution up to t = 10 5 yr. Solid particles with rz. Js 
1 cm, 10 cm, and 10 3 cm are all short-lived and their evclu 
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Fig. 5. Radial distributions of surface density of solids at selected times (t = 2 x 10 4 yr, t = 6 x 10 4 yr, and t = 10 5 yr) for SL particles. 
The left-hand panels correspond to the evolutionary scenario with a = 0.01, whereas the right-hand panels correspond to the scenario with 
a = 0.001. Dotted lines show the surface density of the gas x 10 -2 . Solid lines show the surface density of 1 cm particles, dashed lines show 
the surface density of 10 cm particles, and dash-dotted lines show the surface density of 10 3 cm particles. 


tion is shown on Fig. 5. If solid material is concentrated in such 
particles, it decouples quickly from the gas and is lost to the 
central star in a relatively short period of time. The case of 10 
cm particles offers the extreme example. Such particles are lost 
in less than a few x 10 4 yr regardless of the magnitude of a. 
Solid material concentrated into 10 3 cm particles is lost in about 
10 5 yr, whereas solids concentrated into 1 cm particles are all 
accreted in a few x 10 5 yr. None of these SL particles survives 
in an average sense, until t = 10 6 yr. As can be seen from Fig. 5 
radial distribution of the surface density for SL particles can be 
quite complicated. Consider, for example, 10 3 cm particles dur- 
ing the evolution characterized by a = 0.001. Although most of 
them accrete quickly onto the central star, those initially located 
at the outer edge of the disk linger there for a relatively long time 
before starting their inward movement. This rather intricate be- 


havior is the result of large radial and temporal changes in the 
aerodynamic regime encountered by particles. 

5. Discussion 

This paper represents an attempt to study the global evolution 
of solids in a viscously accreting protoplanetary disk. This is a 
first venture into a complex problem, and we have accordingly 
Simplified the problem to mere essentials. Two central issues for 
modeling global evolution of solids are aerodynamic forces and 
coagulation processes. We have developed a numerical method 
to study the effects of aerodynamic forces. Adding a coagulation 
process to our numerical model will allow us to test whether the 
body of presently embraced theoretical concepts (outlined in 
Sect. 1), when applied all at once, leads to the formation of a 
planetary system that looks anything like our own. This will be 
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presented in the next paper. Despite the obvious limitation ot 
studying aerodynamic forces alone, they yielded a number ol 
interesting conclusions. 

Perhaps the most interesting (and upon reflection the most 
obvious) result of our calculations is that solid particles that are 
larger than 0. 1 cm are not really entrained inthc gas. This finding 
can be best illustrated by following the total mass of solid parti- 
cles in a disk as a function of time. Fig. 6 shows such a total mass 
evolution for particles of different sizes. Examination of Fig. 6 
makes it clear that solids and gas evolve differently. This seems 
to preclude any cosmogonical theory that postulates that the ra- 
dial distribution of solids remains in any particular relationship 
to the momentary radial distribution of the gas. In particular, our 
results suggest that the broadly used concept ot the minimum 
mass solar nebula, in which the radial distribution of the gas 
in the disk is reconstructed from the radial distribution of the 
solid material as found in the present-day solar system, is con- 
ceptually Hawed 5 . As our solar system is readily accessible to 
observation, it may seem quite salutary to use any data we have at 
hand to advance our knowledge about the progenitor solar neb- 
ula as well as other protoplanetary disks. However, our results 
show that such an approach, although intellectually appealing, 
is futile, as it must inevitably lead to spurious results. Instead 
of trying to reconstruct a protoplanetary disk from a planetary 
system, we should simulate the formation of a planetary system 
by means of comprehensive evolutionary calculations. 

Weidcnschilling (1988) discusses a variety of coagulation 
processes that are relevant for solar nebula conditions. Natu- 
rally, if a solar nebula exemplifies contemporaneously observ- 
able disks, these same processes are responsible for coagulation 
in all protoplanetary disks. For particles of the sizes considered 
here, the predominant mechanisms are coagulation due to ra- 
dial motion, vertical settling, and turbulence. For our models 
of an accretion disk, characterized by a relatively small value 
of a, turbulent coagulation is generally slower than the coag- 
ulation due to the radial motion of particles. Vertical settling, 
which gives the fastest coagulation in the quiescent disk, is not 
very effective in the turbulent disk. Altogether, for the purpose 
of the present discussion, we may assume that coagulation is 
dominated by radial drift. The mass of the particle increases 
at the rate dm/dt ss 7rs 2 p3AV;, where AV r is assumed to be 
of the order of \V r d\- Thus, the characteristic coagulation time 
is t c = ml (dm/dt) « (p s /pd)(s/ Al£). Using this formula, 
one can estimate a characteristic time scale of coagulation for 
various sized particles at different times and different disk lo- 
cations. For example, consider a disk at t = 10 5 yr. Our model 
with a = 0.01 gives ~pl ~ 10" 13 g cm" 3 at 10 AU from the 
star for small and large particles, and ~pl ~ 10" 14 g for medium 
particles. Assuming p s of about 1 g cm" 3 and using Figs. 1 - 
3 to estimate the value of AV r , we calculate t c ~ 10 3 yr for 


5 As the name suggests, the minimum mass solar nebula concept was 
originally used to estimate the total minimum mass of the disk, which 
is a perfectly viable exercise. Later, however, this simple concept was 
extended (see for example Hayashi et al. 1985) to construct models of 
protoplanetary disks giving radial distribution of E, H, 7\ and p. 


small particles, t c « 10 4 for medium particles and t c ^ 1C 
yr for large particles (here we used AV$ instead ot AV r to c> 
timate £ c , because \A%\ > |Al£| for large particles). Thes ■ 
estimations underscore the necessity of our global approach t > 
the evolution of solid matter in a protoplanctary disk, especiall 
for medium and large particles. The coagulation time scale fc r 
small particles is indeed short in comparison with their cha: 
acteristic radial evolution time scale. A small particle move , 
only about 0.1 AU from its original location at 10 AU befoi : 
it becomes a part of a larger particle. However, medium pa - 
tides have coagulation time scales comparable to their radi 1 
evolution time scales. An average medium particle may trav. 1 
all the way to the central star during its lifetime. This is also tri e 
for large particles, which take a long time to coagulate. A large 
particle has very slow radial velocity, but it can travel a long 
distance before coagulationg into even bigger particle. Thus, a 
local approach may work for small particles that are slow ar d 
coagulate quite quickly. Medium particles coagulate quickly b t 
are very fast, whereas large particles move slowly but they exi t 
for a long time; in both cases global calculations are needed o 
calculate the radial distribution of solid matter. 

Considering a turbulent instead of a quiescent disk is ir..- 
portant for the physical conditions describing the gas and f r 
aerodynamic forces and coagulation processes. The properti s 
of the turbulent geseous disk changes with time, thus changii g 
aerodynamic regime encountered by particles. Even if coag i- 
lation of small particles happens fast, the subsequent evolutic n 
of large particles takes place on much longer time scales, ov t 
which it is necessary to take into account changing gas conci- 
tions. Our model assumes that turbulence is present throughout 
the entire phase of solids evolution. This point of view co i- 
trasts with the often-expressed opinion that turbulence decayed 
soon after the end of parent cloud collapse (Hayashi et al. 1985; 
Safronov, 1991). Only in thecase when turbulence is exclusively 
associated with convection can one expect that a disk becomes 
laminar after particles with sizes < 10" 2 cm coagulate, mak- 
ing a disk transparent and stable against convection. For this to 
happen the coagulation process must be very efficient, leavi g 
no small particles behind. However, infrared, optical, and near 
ultraviolet observation of T Tauri stars show the presence of ac- 
cretion (see Bertout et al. 1988) for as long as there is eviden e 
for a disk being present around a star. We may thus infer tf at 
protoplanetary disks maintain turbulence throughout their liie- 
times. Such turbulence may be due to convection (if the process 
of small particles coagulation is inefficient), or may be the re- 
sult of magnetorotational instability (Balbus & Hawley, 199F 
or even purely rotational instability (Dubrulle, 1992). There 
fore, in our opinion, a disk model with “generic” turbulen :c 
parameterized by constant a is a more realistic approximati 
of real protoplanetary disk than a laminar disk. 

Besides being important for the evolution of the gas, the i 
elusions of turbulence has a direct consequences for dynamic: 
of solid particles. By virtue of being coupled to the turbulent gas 
particles are subject to “pressure”, “viscosity”, and “diffusivit v.’ 
Obviously, the extent to which these effects are important oe 
pends on the effectiveness of gas-solids coupling. For examp le 
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Fig. 6. Summary of solid particle evolution. The panels show the total mass of particles in the disk as a function of time: the left-hand panel 
corresponds to the evolutionary scenario with a - 0.0 1 , and the right-hand panel corresponds to the scenario with a = 0.00 1 . The correspondence 
between line styles and particle sizes is the same as in Figs. 5 and 6. In addition, lines are labeled by a number indicating particle size in cm. 


the radial velocity of a 0.1 cm particle (which is well coupled 
to the gas) in an a = 0.0 1 turbulent disk is about 2 orders of 
magnitude larger than the radial velocity of the same particle in 
a hypothetical quiescent disk characterized a-disk gas density 
and temperature. Other things being equal, the turbulent disk 
viscosity causes particle radial motion much larger than the ra- 
dial motion of particle in a quiescent disk, which is exclusively 
due to the drag resulting from existence of the relative tangen- 
tial velocity between a particle and the gas. Medium and large 
particles are less coupled to the gas and their radial velocities 
are dominated by the drag and not by turbulent viscosity. 

Inspecting Fig. 6 we conclude that the total mass of 10 — 10 3 
cm particles drops below the total mass of solid material in the 
present-day solar system (~ \0~ 4 Mq) on a time scale of about 
IO 4 - t0 5 yr. Therefore, the coagulation time scale for particles 
of such sizes must be faster than 10 4 yr, otherwise there will be 
no material left to form the planets. Above we have estimated 
that t c ~ 10 4 yr for 1 cm particles located at 10 AU after 10 5 
yr from the beginning of our calculations. Extending this esti- 
mation for all 1 cm particles, we may conclude that they will 
increase their size and slow down before spiraling into the cen- 
tral star. However, the value of t c is sensitive to all of the above 
assumptions, and the question of the solid matter survivability 
cannot be answered without including the coagulation process 
into our calculations. 
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